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ABSTRACT 

d. 

Aims. A pilot project has been carried out to measure circumstellar CO emission from three OH/IR stars close to the Galactic centre. 
(~ l , The intention was to find out whether it would be possible to conduct a large-scale survey for mass-loss rates using, for example, the 

^ |. Atacama Large Millimeter Array (ALMA). Such a survey would increase our understanding of the evolution of the Galactic bulge. 

I ' Methods. Two millimetre- wave instruments were used: the Nobeyama Millimeter Array at 115 GHz and the Submillimeter Array 

P ' at 230 GHz. An interferometer is necessary as a 'spatial filter' in this region of space because of the confusion with /n/erstellar CO 

emission. 



^ . Results. Towards two of the stars, CO emission was detected with positions and radial velocities coinciding within the statistical 

^ I ' errors with the corresponding data of the associated OH sources. However, for one of the stars the line profile is not what one expects 

for an unresolved expanding circumstellar envelope. We believe that this CO envelope is partially resolved and that this star therefore 
, is a foreground star not belonging to the bulge. 

^ ■ Conclusions. The results of the observations have shown that it is possible to detect line profiles of circumstellar CO from late-type 

' stars both within and in the direction of the Galactic bulge. ALMA will be able to detect CO emission in short integrations with 

' sensitivity sufficient to estimate mass-loss rates from a large number of such stars. 

' Key words. Stars: AGB and post-AGB; (Stars:) circumstellar matter; Stars: mass-loss; Galaxy: bulge; Radio lines: ISM; Techniques: 

interferometric 

(N 
O 

0^ 1 . Introduction several OH/IR stars close to the GC but detected CO emission 

O ■ from only one star: OHO.3-0.2 (Baud et al., 1975). This star 

- ■ Rapid stellar mass loss occurs at the so-called red giant and, j^^^ ^ ^-^ ^^^-^^ ^^j^^j (_34^ j^^-i^ therefore the 

^ in particular, the asymptotic giant branches (RGB and AGB). circumstellar emission is not affected by the interstellar back- 
Due to high stellai- density and rapid star formation, the Galactic .^^^^^ g^^j^^i^j^ ^j^^^ ^^^^^^ i^^jj ^j^^ ^^j^^^ _2oo to 
bulge contains large numbers of RGB and AGB stars. It is im- ^200 km s ' . None of the other candidate stars were detected 
portant to measure the mass-loss rates of these stars for a picture ^^^^^^^ confusion with interstellar emission, 
of the recirculation of matter and metal enrichment in this re- 
gion of the Milky Way. In addition, stars in the central bulge are The present project employs a different observing technique 
at distances that differ from one another by merely a few percent in an attempt to detect the circwmstellar CO emission in the midst 
leading to an accurate estimate of the mass-loss-rate distribution, of interstellar CO emission. A radio interferometer with suitable 
The most accurate method of estimating the stel- baselines can be used as a 'spatial filter' by resolving most of the 
lar mass-loss rates is based on the CO rotational spec- interstellar background but leaving the circumstellar emission as 
tral lines (cf. Ramstedt et al., |2008] l. The first attempt unresolved point sources. 

to detect such cii-cumstellar lines in the vicinity of the ^o determine optimal baseline lengths and the most 

Galactic centre (GC) was made by Mauersberger et al. favourable CO lines, we started a pilot experiment. We chose 

CMi using the IRAM 30-m telescope. They detected the ^^^^^ q^/jr ^^^^^ ^j^^g ^j^^ p^^iji^^ with strong 

7=1^0 and 2^1 CO hnes from the proto-planetary neb- jp, fl^^^g^ ^j^^ low-to-moderate radial velocities. In 2003 

ula (PPN) OHO.9+1.3 which has a high radial velocity _ 2004 we used the Nobeyama Millimeter Amy (NMA) at 

(-1 10 km s-i ) making it likely to be physically close to the (^O, 7= 1 ^0) and in 2005 we used the SubMillimeter 

GC. They also proposed to look for CO emission from OH/IR ^^-ay (SMA; Mauna Kea, Hawaii) at 230 GHz (7=2^1). 
stars, since such stars are believed to be the progenitors of PPNe. 

Winnberg et al. (1991 ) used the same radio telescope to observe This research note presents the main results of these two data 

sets and outlines the prospects for future systematic surveys of 

Send offprint requests to: A. Winnberg late-type stars in the Galactic bulge using the Atacama Large 
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Millimeter Array (ALMA). A conference report of this project 
appears in Winnberg et al. ( 120061 1. 
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Fig. 1. Circumstellar CO(2-»l) line profile of OHO.3-0.2. The 
thick 'staircase' line is the observed spectrum; the thick curved 
line is a fitted parabola; the dotted line depicts the residuals (i.e., 
observed spectrum minus parabola). The arrows indicate the ve- 
locities of OH masers. 
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Fig. 2. Circumstellar CO(2— >1) line profiles of 
OH359. 762 +0.120. The thick 'staircase' line is the observed 
spectrum with baseline lengths shorter than 25 k/1 excluded. 
The parabola, the residuals and the OH maser velocities are 
represented as in Fig. [T] In addition, a thin 'staircase' line 
indicates an observed spectrum with all projected baseline 
lengths included. 



2. Observations 

The first observations in this project were done with the 6- 
element array (NMA) at Nobeyama Radio Observatory, Japan, in 
November 2003 and in January 2004 at 1 15 GHz. Three OH/IR 
stars were selected for observations on the basis of strong IR 
fluxes: OH359. 1 17 -0.169 , OH359.762+0.120 and OH359.971- 
0.119 (Ortiz et al., 120021 ). The relevant data for these stars are 
listed in Table[T] where f [15] denotes the flux density at a wave- 
length of 15 yum as measured by the Infrared Space Observatory 



(ISO). A clear signal was obtained from OH359.97 1-0.1 19 (Fig. 
[3] dotted line). The position of the CO line agreed with the po- 
sition of the 1612-MHz OH line. However, the line width was 
too narrow for an unresolved circumstellar CO source and the 
radial velocity of the line was close to that of the red-shifted OH 
line component rather than to the average velocity of the two OH 
components. 

Data were taken using two digital correlation spectrometers 
simultaneously. One of them consisted of 128 channels with a 
width of 8.0 MHz. The other one had 1024 channels of width 
31.25 kHz. 

In June 2005 the same three stars plus the previously de- 
tected star OHO.3-0.2 were observed using the 8-element array 
(SMA) on Mauna Kea, Hawaii, at 230 GHz. The spectrome- 
ter was configured to give a resolution of 3.25 MHz. The test 
star OHO.3-0.2 was detected and showed properties in accor- 
dance with the single-dish data (Fig. [T]). In addition the stars 
OH359.762-h0.120 and OH359.971-0.119 were detected (Figs 
|2]and[3l). 

Typical synthesized beam FWHM values of images pro- 
duced with all basehnes were 4" x 7" for the NMA and 3" x 
4" for die SMA. 



3. Data reduction 

3.1. The NMA 

The data were calibrated using the Nobeyama internal pro- 
gramme package and subsequently the Astronomical Image 
Processing System (AIPS) was used in a standard way to obtain 
images and spectra. 

We found strong ripples in the images that could be elim- 
inated by removing data from projected baselines shorter than 
10 k/1 (26 m). Due to atmospheric phase instability, we decided 
to remove all projected baselines longer than 40 k/l (104 m) and 
to introduce a gaussian baseline-length taper such that baselines 
of length 40 k/l got a weight of 30 %. 

Maps were made with 256 x 256 pixels of size 0.5" with 
uniform weighting. They were 'cleaned' using the standard 
Hogbom/Clark algorithm with a gain of 0. 1 and a minimum flux 
density per clean component being the product of the beam dy- 
namic range (l/|strongestsidelobe|) and the expected rms noise 
fluctuations. This ensured that no 'overcleaning' took place in 
the rather noisy maps. 

Data cubes were made consisting of images from the 20 
central channels of the low-resolution spectrometer covering 
416 kms For the high-resolution spectrometer similar data 
cubes were made consisting of 256 channels where each chan- 
nel resulted from averaging 4 original channels. Therefore these 
spectra covered 83.2 km s"' with a resolution of 0.325 km s"' . 

3.2. The SMA 

Data from this instrument were retrieved using software from 
the Radio Telescope Data Center (RTDC) of the Center for 
Astrophysics (CfA). Calibration was performed using the im- 
age processing package IDL-MIR at Academia Sinica Institute 
of Astronomy and Astrophysics (ASIAA). Further data reduc- 
tion was made using both AIPS at Onsala Space Observatory 
(OSO) and Miriad (SMA version) at ASIAA. 

The uv data were investigated for the presence of spa- 
tially extended emission by plotting the visibility amplitude 
as a function of projected baseline length. Based on such 
plots it was decided to exclude baselines shorter than 25 k/l 
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Name 


RA 


Dec 




* exp 


^[15] 


Ref. 




(J2000) 


(J2000) 


(km s-') 


(kms-') 


(Jy) 




OHO.3-0.2 


17'' 47"' 06.95'' 


-28° 44' 42.2" 


-341.0 


14.5 




1 


OH359. 117-0. 169 


17'' 47"" 21.79" 


-29° 47' 42.2" 


-88.5 


21.2 


6.9 


2 


OH359.762+0.120 


17" 44"" 34.95'^ 


-29° 04' 35.2" 


-5.7 


15.3 


13.8 


3 


OH359.971-0.119 


17'' 46"" 00.94" 


-29° 01' 23.6" 


-8.5 


19.3 


9.5 
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References: l.-Fix & Mutel ( fT984t ; 2.-Sevenster et al. ( fTWTl i; 3.-Lindqvist et al. ([1992) 
Table 2. Results 



Name 


Trans. 




RA 


Dec 








M 




(J) 




(J2000) 


(J2000) 


(kms-') 


(kms-') 


(Jy) 


(Moyr-') 


OHO.3-0.2 


2-^1 


17'' 


47'" 06.978^(0.009) 


-28° 44' 42.8" (0.2) 


-340.1(0.6) 


14.0(0.8) 


1.8(0.1) 


5 X lO-'* 


OH359. 117-0. 169 


1^0 












<0.2 






2^1 












<0.54 




OH359.762+0.120 


1^0 












<0.15 






2^1 


17'' 


44'"34.979''(0.012) 


-29° 04' 36.6" (0.2) 


-11(1) 


13(1) 


0.6(0.1) 


4 X lO-'' 


OH359.971-0.119 


l->0 


17'' 


46'" 00.75" (0.04) 


-29°01'21.5"(0.7) 






-0.4 






2^1 


17'' 


46'" 00.909" (0.008) 


-29° 01' 23. 2" (0.2) 
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(32.5 m) for OH359.762+0.120 and shorter than 30 kA (39 m) 
for OH359.971-0.1 19 in order to avoid, as far as reason- 
able, contamination by residuals of interstellar emission. For 
OH359. 117-0. 169 no evidence of significant interstellar emis- 
sion was found. No interstellar emission was found in the IF 
band of the test source OHO.3-0.2 as expected. 

Map-cubes were made with 256 x 256 pixels of size 0.5", 
however, this time with natural weighting. They were 'cleaned' 
using the standard Hogbom/Clark algorithm in a manner similar 
to the treatment of the NMA data. When a compact source was 
seen in one of the channels within the OH velocity span at a 
position close to the position of the OH source, we fitted a two- 
dimensional elliptical Gaussian and extracted a spectrum at the 
pixel closest to this position. 

After exclusion of the shorter baseline lengths all the maps 
are free from interstellar CO emission except for occasional 
point sources. We interpret them as unresolved rests of inter- 
stellar emission. Since most of the emission is resolved there 
are 'negative point sources' as well, i.e. unresolved 'dips' in the 
interstellar background. The positions of the OH/IR stars are 
known to 1" or better, and therefore there is never any doubt 
about the identification of a circumstellar CO source. In addition 
to the position evidence, there is the radial velocity evidence that 
strengthens the identification case quite considerably. 



4. Results 

Table |2] lists the measured parameters of the CO sources as- 
sociated with the four stars observed. The 1-cr errors are 
given in parentheses after the values. For OHO.3-0.2 and 
OH359.762-h0.120 (2^1) least-square fits of parabolas to the 
CO line profiles have been made, assuming that the data are from 
unresolved, optically thick emission; 
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Fig.3. Circumstellar CO(2^1) line profiles of OH359.971- 
0. 1 19. The definitions of the lines are given in the caption of Fig. 
12] except in this figure the dotted 'staircase' line is the CO(l— >0) 
spectrum as observed with the NMA. 



where 5 m is the maximum flux density at the radial velocity 
Vm, and Ve is half the line width at zero intensity. These three 
parameters are given in the table together with the statistical 
errors from the fitting procedure. No CO emission associated 
with OH359.1 17-0.169 (l-»0 and 2^1) and OH359.762+0.120 
( 1 — >0) was found and 3-cr upper limits are given for Sm- The CO 
line profiles for OH359.971-0.119 (1^0 and 2^1) are incom- 
patible with a parabola, although the CO positions are coincident 
with the OH position, and only approximate values of the maxi- 
mum flux densities are given (Fig.[3]l. 

OH359. 117-0.169 was not detected at any of the two CO 
lines and we do not know the reason for it. One guess would be 
that this is due to heavy absorption by interstellar CO clouds in 
front of the star. Such a case could possibly occur for a star that 
is situated at a distance beyond the GC. 
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OH359.762 +0.120 was detected with a rather poor S/N. 
Therefore there are quite large errors associated with the ele- 
ments of the fitted parabola (Fig. |2] and Table |2l). Within these 
errors, the mean radial velocity of the line and the line width are 
compatible with the systemic velocity of the star and its enve- 
lope expansion velocity as measured from the OH line profile. 
Notice that the exclusion of short baselines did not improve the 
detection of the line significantly. It merely improved the spec- 
tral baseline. 

OH359.97 1-0.1 19 was detected with moderate to good S/N. 
However, the CO line is narrow and close to the 'red-shifted' 
OH line component, i.e. the backside of the envelope (Fig. |3]l. 
A similar, but considerably weaker, line was observed with the 
NMA at 115 GHz (dotted line). Notice that the line profile ob- 
tained when all the uv data were included (thin solid line) shows 
a much stronger and broader 'red-shifted' line and even a 'blue- 
shifted' counterpart. The 'negative signal' near -20 kms"' is 
probably caused by a 'dip' in the interstellar background, as dis- 
cussed above, but in this case it is resolved. 

Mass-loss rates have been calculated for OHO. 3-0.2 and 
OH359. 762 +0.120 (2— >1) using a new formula based on the 
original equation by Knapp & Morris ( IT985) but containing 
constants determined by least-squares fits to physical models 
(Ramstedt et al., I2008I I. and the values have been entered in 
Table |2] For these calculations a distance of 8 kpc to the GC 
(Reid [T993] l and a CO /H2 abundance ratio of 2 x 1 (Ramstedt 
et al.. l2008l l have been assumed for both stars. The values of Ve 
have been taken as the expansion velocities of the CSEs. Both 
mass-loss rates are normal for OH/IR stars (we caution that the 
validity of the Ramstedt et al. formula has only been tested in 
the mass-loss-rate range 10"^ to 10"^ Mq yr"'). 

5. Discussion 

5.1. OH359. 762+0. 120 

There is little doubt, in spite of the poor S/N, that the detected 
CO source is a true circumstellar source associated with the 
OH/IR star Because of the high bolometric magnitude of this 
star and its strong OH emission, some people doubt the asso- 
ciation of this star with the Galactic bulge (see discussion by 
Blommaert et al.. |1998l l. However, as pointed out by Blommaert 
et al. ( 119981 1. the OH radiation is heavily scattered by interstel- 
lar free electrons (Frail et al., 119941 1 making it very likely that it 
resides close to the GC. Our result of a weak CO source lead- 
ing to a normal mass-loss rate for an assumed distance of 8 kpc 
supports this conclusion. 

5.2. OH359.971-0.119 

This star (as well as OH359.762+0.120) has detected 43-GHz 
SiO emission (Lindqvistet al.. 119911 ). However, the v-2, 7=1— >0 
SiO line is much stronger than its v= 1 counterpart. This suggests 
a cool dust temperature and a high mass-loss rate (Nakashima & 
Deguchi,|2007ll. 

The CO source associated with this star, on the other hand, 
is an enigma. The line profile is similar to that expected from a 
well-resolved circumstellar envelope where unresolved emission 
is left at the front and back sides (see for example the central CO 
source in U Cam, Lindqvist et al. . 119991 1. For such a picture to be 
true, the star needs to be at a small distance and would not belong 
to the Galactic bulge. For example, assuming that the diameter of 
the CO envelope is 2 x 10'^ cm (which might be an overestimate 
for the 2— >1 transition) and that the angular diameter is 10", the 



distance would be only about 1.4 kpc. The IR properties of this 
star also are such that it is arguable whether it belongs to the 
bulge (cf. Ortiz et al.. l2002b . 

Another explanation of the line profile would be that it is 
heavily distorted through strong absorption of interstellar CO 
in front of the source. A third possibility - although improba- 
ble - is that the physical conditions in this CSE are such that 
they favour weak maser action along radial directions (Morris, 
1980). Finally, there remains the (improbable) explanation that 
the source is an unresolved remnant of interstellar CO emission 
that happens to be at the same position and radial velocity as 
OH359.971-0.119. 

Observations of higher-energy CO lines might be a possible 
way of revealing the true nature of this source. However, given 
the available data, we favour the first alternative, i.e. that we have 
found a relatively nearby OH/IR star whose CO envelope is re- 
solved by both arrays used. The salient points that support this 
conclusion are: 

- There are two CO line components whose positions and ra- 
dial velocities correspond to those of the OH line compo- 
nents. 

- The redshifted CO line component grows stronger and the 
blueshifted line component emerges when all baselines (in- 
cluding the short ones) are included in the imaging. 

- There is a weak line component at 115 GHz whose radial 
velocity coincides with that of the redshifted 230-GHz com- 
ponent and whose position coincides with that of the OH/IR 
star. 

- The 1 15-GHz line is much weaker than the 230-GHz coun- 
terpart which mostly is the case for circumstellar CO lines in 
OH/IR stars especially those with cool CSEs (Heske et al., 
fT990l ). 

6. Conclusions 

Our pilot project has shown that it is possible to detect circum- 
stellar CO envelopes of OH/IR stars close to the GC: out of 
three stars selected, one was detected at both 115 and 230 GHz 
(OH359.971-0.119) and another one was detected at 230 GHz 
only (OH359.762+0. 120). OH359.971-0. 1 19 probably does not 
belong to the Galactic bulge, but this fact is irrelevant for the is- 
sue at stake - this star too would have been hard to detect using 
a single-dish telescope. 

We would have liked to observe the same three stars at 
345 GHz (7=3— >2) and to try out somewhat longer baselines, 
but both these requests require excellent atmospheric conditions 
and therefore the competition for observing time is strong. 

We have no doubt that ALMA will be able to detect a large 
number of OH/IR stars in the inner bulge (Olofsson. 120081 ) . 
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